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We show that the dissipationless spin current in the ground state of the Rashba model gives rise 
to a reactive coupling between the spin and charge propagation, which is formally identical to the 
coupling between the electric and the magnetic fields in the 2 + 1 dimensional Maxwell equation. 
This analogy leads to a remarkable prediction that a density packet can spontaneously split into 
two counter propagation packets, each carrying the opposite spins. In a certain parameter regime, 
the coupled spin and charge wave propagates like a transverse "photon" . We propose both optical 
and purely electronic experiments to detect this effect. 
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The generation and manipulation of spin current is es- 
sential to the rapidly developing field of spintronics^. For 
this purpose, the analogy between photonics and spin- 
tronics greatly helped our conceptual developments, and 
this analogy lead to the celebrated Das-Datta proposal 
for a spin-field transistor^. However, earlier attempts to 
realize this conceptual device were based on the rather 
incomplete analogy between the electron spin and the 
photon, and were plagued by many issues such as low 
spin injection rate and the requirement of the ballistic 
spin transport. 

Recently, the remarkable phenomenon of the dissipa- 
tionless spin current has been theoretically predicted^ . A 
electric field generates a spin current described by the 
response equation 

jj = asEi-jkEk (1) 

where is the current of the i-th component of the spin 
along the direction j, Cijk is the totally antisymmetric 
tensor in three dimensions and the spin Hall conductiv- 
ity (Ts does not depend on impurities. Since both the spin 
current and the electric field are even under the time re- 
versal, this equation describes a reactive response which 
does not dissipate energy. One natural consequence of 
this equation is the intrinsic spin Hall effecti^-, which has 
been recently observed experimentally in the hole doped 
systems^. Another consequence is the dissipationless spin 
current in the ground stated. In the above equation Q), 
the electric field can be either externally applied, or can 
be spontaneously generated in systems without inversion 
symmetry. In a two dimensional electron gas (2DEG), 
the confining potential along the z direction breaks the 
inversion symmetry, and leads to a internal electric field 
Ez in the ground state. According to equation (QJ, there 
is a spin current in the ground state, jj = jo^ij, where 
eij is the antisymmetric symbol in two dimensions with 
hj = x,y. 

In this work, we shall show that the dissipationless spin 
current in the ground state makes the analogy between 
photonics and spintronics formally exact. In 2 -I- 1 di- 
mensions, the electric field has two components, while 
the magnetic field has only one component. If one iden- 
tifies them with the in-plane components of the spin den- 



sity and the charge density, respectively, the Boltzmann 
transport equation for the coupled spin and charge wave 
is formally the same as the Maxwell equation describ- 
ing the electromagnetic fields, where the "speed of light" 
is given by the Rashba coupling constant. This behav- 
ior is in sharp contrast to the conventional Boltzmann 
equation for the decoupled spin and charge dynamics in 
semiconductors, where only purely diffusive, but no prop- 
agating motion is predicted^ . The photonic analogy helps 
our understanding on how density gradient and time de- 
pendence can generate spin density, and leads to many 
novel predictions. We shall show that there is a parame- 
ter regime, reachable experimentally, where the coupled 
spin-charge wave propagates as a under-damped "pho- 
tonic mode". A density packet will split spontaneously 
into two counter propagating packets, each carrying the 
opposite spins. This mechanism enables injection of spins 
and spin currents. The Boltzmann transport equations 
for the Rashba model have been studied previously in 
the diffusive regioniSili. The coupled spin-charge wave 
propagation is a new result of this work. 

A spin 1/2 Hamiltonian which includes spin orbit cou- 
pling can be written in the following general form: 

H=^ + X'{p)a\ i = x,y,z (2) 
2m 

where A* (p) is a odd function of p, in order to preserve 
the time reversal symmetry. This includes a wide range 
of spin-orbit couplings, including 2D Rashba and Dres- 
selhaus couplings and the 3D spin splitting of the conduc- 
tion band in strained semiconductors^. The phase space 
density distribution function n p {p, r, t) and the energy 
matrix tp{p^r^t) are 2x2 matrices, and can be decom- 
posed as: 

npip, r, t) = r, t) + S\p, r, t)a^ 

eF{p,r,t) ^ es{p,r,t) + el{p,r,t)a\ i = x,y,z (3) 

In this letter we consider the system in the absence of 
external fields, such that e_F(p, r, i) = p^ /2m + X^{p)a^. 
The influence of electric and magnetic fields on the sys- 
tem is described in a future longer publication^. The 
Boltzmann equation reads: 
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where the right hand side is the colhsion term expressed 
in the relaxation time approximation, and n'p is the equi- 
hbrium value of nF(j>Tr,t). t is the momentum relax- 
ation time. Although this approximation does not take 
into account the self-energy effects, it turns out to be 
qualitatively and quantitatively correct, as we see from 
comparison with the solution involving the self energy in 
some special cases of spin-orbit couplingi2iii. We trace 
out the matrix dependence of the distribution function as 



well as that of the energy, and integrate the continuity 
and the current equations over the Fermi volume^ . After 
linearization, we obtain^: 

dnjr, t) _^ dMr, t) dX' (p) dS' (r, t) 
dt drf dpi dn 



dt 



r,t) _ ^d'^S^{r,t) d\^{p)dn{r,t) 



dpt dn 



4mD dX'dS^{r,t) 



'dpr 



I 



2m 
IT 



D 



dpr dpr ' 



dpr dpr 



S%r,t) 



(6) 



where D 



2m2 



is the diffusion constant and fi 



is the mobility, and the carriers are electrons of charge 
— e. Aside from the self energy renormalizations, Eq.jS] 
gives the same result asiSiii when particularized to the 
Rashba-spin-orbit coupling. The last term in the spin 
continuity equation represents the spin relaxation due 
to Dyakonov-Perel (DP) mechanism^^. In the case of 
Rashba systems, the spin orbit coupling is A* = aeijzPj, 
a has the dimension of velocity, and the continuity equa- 
tions become (from now on we use di — d/drt and 
dt = d/dt) : 

dtn = Dd'in-ae^'''d^S^ 

dtS'' = Dd^S'' - ae'^'^d^n + ^j^{5kzd^S' - dkS') - 



(7) 



where — {x,y) {i — 1,2) since charge and spin mo- 
tion is now confined entirely to the 2D plane. Within 
the current microscopic approximation, the DP spin re- 
laxation time is given by = (^^)^I?, however, in 
the subsequent discussions, we shall treat D and as 
independent, phenomenological parameters. 

Let = {S'',S^), 11 = x,y,z, k = x, y. We can write 
the two dimensional vector S^{r,t), k = x,y in the most 
general form as a sum of a longitudinal vector 5^ (r, t) 
and a transversal vector S'^(r, t): 

5'= = 5i + 5^; dkSl^ = 0, e^jd^Si = (8) 
Substituting this decomposition into Eq. Q , we find two 



sets of coupled equations: 

dtn = Ddfn - aekidiS^ 
dtSl^ = Dd^S^-aeMd,n^ 
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Hence the charge density couples only to the transverse 
spin component, while couples only to the longitudinal 
spin component in a purely diffusive fashionii. In the 
spin continuity equation of Eq. ((SJ, we see that the a 
term is nothing but the divergence of the dissipationless 
spin current in the ground state = aekiU. We shall 
see that this term plays the crucial role leading to the 
coupled spin-charge propagation. 

At this point we come to a remarkable realization that 
the Boltzmann equation (PJ for the coupled charge and 
transverse spin transport is exactly the Maxwell's equa- 
tion in 2 + 1 dimensions! In order to facilitate the com- 
parison, let us first focus on the large spin-orbit coupling 
limit, where we neglect the D and the terms in Eq. 
0. In 2 -I- 1 dimensions, the source- free Maxwell equa- 
tions are given by 



e'^-'Pd^F.p = 



(11) 
(12) 
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where = (0,1,2/). In 2 + 1 dimensions, the magnetic 
field has only one component, given by — F^y, and the 
two components of the electric field are given by Ei = Fqi. 



If we make the identification : 



Ei and a 



c, we see that the three Boltzmann transport equations in 
Eq. are exactly the three Maxwell's equations in the 
vacuum of 2 + 1 dimensions, namely the Faraday's law 
of induction, the Ampere-Maxwell law, and the Gauss's 
law. More generally, the -^S!^ term can be interpreted as 
light propagation in a metallic media, with the current 
density in the Ampere's law given by the Ohm's law, and 
the D terms can be interpreted as due to light diffusion 
in a random media. 





We see that an initial density wave packet spontaneously 
splits into two counter-propagating packets, each carry- 
ing the opposite spin. This phenomenon can be elegantly 
interpreted in the "photonic" language. In 2 -f 1 dimen- 
sions, the magnetic field is always pointing along the ~\-z 
direction. Since the propagation vector k, being propor- 
tional to the Poynting vector, is given by k oc E x B, it 
uniquely determines the direction of the transverse elec- 
tric field. Translating from the "photonic" language into 
the "spintronic" language, we see that the mode prop- 
agating along the -l-a; has spins along the +y direction, 
while the mode propagating along the —x has spins along 
the —y direction. The split wave packets carries a spin 
current JJf, which is a reflection of the spin current in 
the ground state of the Rashba model. For a simple esti- 
mate, a = 3 X lO^m/s, and hence the mode will cross a 
sample of l/im length in 30ps. Considering that the spin 
coherence time in these samples can be larger than Ins, 
it means that the propagation time over l/im distance is 
well shorter than 30-th part of the spin relaxation time 
and can hence be very useful for spin manipulation. 

We now consider the more general situation includ- 
ing diffusion and relaxation. We suppose that a one di- 
mensional stripe of charge density has been created, say 
by transient gratingi^. The initial density is given by 
n{x,y;t = 0) = 5{x). The solution to the full equations 
give S^{x,t) — Sy{x,t) — while S^^Xjt) is generated 
by the spin-orbit coupling: 
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(14) 



{2ttY —{uj — uji){u) — UJ2) 



sy{x,t)^ I I ^-^^__!^i__^e'("*-«^)d^d(z 

(15) 

where uji , LjJ2 are the characteristic frequencies of the sys- 
tem: 



(27r)^ {co — u>i){u) — W2) 



FIG. 1: Charge and spin density for Ts = Ins, a = 3x 10 m/s 
and D = 10~^m^ / s. We see propagation over distances of 
more than lOO/im. Inset: Charge and spin density in the 
diffusive regime, for small values of a ( = Ins, a — lO^m/s 
and D — 10^'^m^/s) has the typical Gaussian decay. 



In conventional theories without spin-orbit coupling, 
the electron transport semiconductors is purely diffusive. 
However, we see that in the limit of strong spin-orbit 
coupling, there is a regime where a propagating, coupled 
spin-charge wave mode is possible. If we neglect the dif- 
fusion and the lifetime terms for the time being, we find 
that the most general solution to the initial condition of 
n{x, y,t ^ 0) = f{x) and S'lj^{x, y, i = 0) = is given by: 

n{x,t) = ]^{f{x + at) + f{x-at)) 

Sy^{x,t) = ^{f(x + at)-f{x-at)) (13) 



wi,2 = i(£'g^ + -^) ± 
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We recognize the propagating mode inside the square 
root. For momenta q > l/2Tsa both characteristic fre- 
quencies contain real parts and hence describe propagat- 
ing waves. However, q must not be as large as to cause 
damping due to the term Dq^t. The condition for this 
gaussian damping to be small is Dq^Ts < 1 for q ~ l/r^a. 
Therefore, the condition for the regime where a propa- 
gating mode could exist is then given by: 



(17) 



This condition can be satisfied in samples where a = 
3 X 10'*to/s and D — IQ^^m? / with longer than 

Insil. In this case, \ — = 10^, much smaller than a. 
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FIG. 2: A modified version of the classic Haynes-Shockley 
experiment. A density packet injected by the emitter spon- 
taneously splits into two counter propagating packets with 
opposite spin. Unlike the settings of the Haynes-Shockley ex- 
periment, one of the two packets propagates to the collector 
without experiencing a sweeping electric field. The time de- 
lay between the injection pulse the the collecting pulse gives 
a purely electric determination of the Rashba spin-orbit cou- 
pling constant. 



In the limit of very long Ts c>o the integrals can be 
solved exactly. We give the expression of in this limit: 



(18) 



The propagating mode is ai ± x = where either one of 
the damping gaussian exponentials becomes unity. The 
spin symmetry is odd in x, the spins propagating in the 
positive and negative x axis directions having opposite 
polarization. Note that for diminishing spin-orbit cou- 
pling a ^ the spin density also vanishes, as it should. 
For finite Ts in a stationary phase-type approximation 
the spin-density solution above gets multiplied by an ex- 
ponential factor exp(— t/2Ts). Impressively, both spin 



and charge can propagate over distances well in excess 
of 100/im and for times well in excess of IOt^, (the full 
time scale is not plotted in Fig^). 

We now propose several experiments to test the cou- 
pled spin-charge wave predicted in this work. One could 
inject the density packet optically, and detect the split- 
ting of the density packet and the associated spin orien- 
tation by optical Kerr rotation. One could also detect 
the spin orientation through the circularly polarized lu- 
minance from the recombination with the majority car- 
riers. Alternatively, one could detect the propagation of 
the density packet purely electrically, by a modified ver- 
sion of the classic Haynes-Shockley experiment'^ . Fig. 
^ describes a narrow sample with light p-doping. Two 
rectifying metal-to-semiconductor point contacts are for- 
ward and reverse biased, respectively, to serve as emitter 
and collector electrodes. After turning on the emitter 
pulse, a electron density packet is injected into the sam- 
ple. In conventional Haynes-Shockley setup, the electron 
packet would be swept to the collector electrode by a 
electric field. In our case, no sweeping electric field is ap- 
plied, but the density packet will spontaneously split into 
two counter propagating packets with the opposite spin 
orientation, with a velocity directly given by the Rashba 
coupling constant a. When the right moving packet is 
captured by the collector electrode, a voltage pulse is reg- 
istered. From the time-delay and the shape of the voltage 
pulse, one can determine the Rashba coupling constant 
and the diffusion constant by purely electric means. This 
experiment illustrates the fact that the injected density 
pulse can take advantage of the spin current in the ground 
state, and propagate without any applied voltage. 
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